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We present the results of analytical calculations of the 3-loop contributions to 
the asymptotic photon vacuum polarization function, in the on shell scheme, 
and of the 4-loop contributions to the on shell QED /3-function. These are 
used to evaluate various 4-loop and 5-loop contributions to the muon anomaly. 
Our analytical contributions to (g — 2) M differ significantly from previous nu- 
merical results. A very recent numerical re-evaluation of 4-loop muon-anomaly 
contributions has yielded results much closer to ours. 
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1. There are several reasons for the increase of interest in analytical multi- 
loop renormalization group calculations, observed during the past decade. Firstly, the 
widespread use of the computer symbolic manipulation systems SCHOONSCHIP [|J, 
REDUCE 0, and FORM |3j stimulates analytical calculations in high-energy physics. 
Secondly, the appearance of elegant methods Q for the evaluation of massless three- 
loop diagrams within dimensional regularization, and their generalization to the massive 
case 0, have made higher-order calculations possible. Thirdly, the detailed study of theo- 
retically interesting and experimentally important QED and QCD predictions entails the 
consideration of the effects of high-order radiative corrections. 

Consequently, several important results have been obtained in gauge theories, for 
example: the analytical 4-loop corrections to the QED /5-function in the modified min- 
imal subtraction (MS) scheme, and in the momentum (MOM) scheme (namely, the \l/- 
function) ||; the 4-loop corrections to cr tot (e + e~ — > hadrons) in QCD 0; the 3-loop 
QED corrections to the photon vacuum polarization function in the MS scheme H; 
the 2-loop corrections to the relation between the on shell (OS) and MS scheme quark 
masses in QCD ||; the 2-loop contributions to OS renormalization constants in QCD and 
QED |10| , |TTJ; the finite parts of 3-loop OS charge renormalization ||. 

These results have been obtained with different computer programs, written in dif- 
ferent symbolic manipulation languages. Specifically, the results of Refs. @-|§ have been 
obtained with the help of the SCHOONSCHIP program MINCER [0, the calculations 



of Refs. 0, [|, [K| were made with the help of several REDUCE packages, whilst the recent 



analysis of Ref. |Tl[ was carried out with the help of the FORM program SHELL2 pj 
Therefore, in view of the importance of the results obtained, it is necessary not only to 
study carefully the final results, but to carry out independent cross-checks of the systems 
and the computer programs involved in the above-mentioned investigations. 

In this work we combine various methods of computer symbolic manipulation Jj]]- 
0, and of evaluating multiloop Feynman integrals 0, |5|], to calculate analytically the 
3-loop QED expression for the photon vacuum polarization function at large q 2 in the OS 
scheme, the 4-loop corrections to the QED /3-function in the OS scheme, and certain 4-loop 
and 5-loop contributions to the muon anomalous magnetic moment a M = (g/2 — 1) M . 
A comparison of our analytical results with recent numerical calculations of a M [|TJ]] is 
presented. 

2. We will start by deriving the basic renormalization group relations between 
QED /3-functions and photon vacuum polarization functions in different renormalization 
schemes. This problem has been previously analysed in Refs. || [15| at the 3-loop level, 



and in Refs. || [T(| at the 4-loop level. Here we will apply the methods of Ref. to the 
4-loop level. 

Let us introduce some basic definitions. Firstly, we define the QED /5-function in 
the MS scheme as 

ain5 m = p.ffi, a) 
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where the renormalized MS scheme coupling constant a is related to the bare one a^ 
by a(fj) = Z 3 (a)aB, \i is the mass unit introduced within dimensional regularization in 
D = 4 — 2e space-time dimensions, and Z 3 (a) is the photon wavefunction renormalization 
constant in the MS scheme. The similar expression for the /^-function in the OS scheme 



reads 



d In a 
d lnm 



-fixed 
-0 



P(a) 



s*er 



(2) 



where a = Z 3 as is the OS scheme coupling constant (the physical fine structure constant), 
m is the electron pole mass, and Z 3 is the photon wavefunction renormalization constant 
in the OS scheme^. 

Using Eqs. (1) and (2), one can obtain the following relation: 
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It can be shown that the ratio a(/u, 
contribution || [15]], i.e. 
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Let us now, following the lines of Refs. [16, 17H , expand the coefficients /3 n , P nJ and g nj in 
powers that correspond to the numbers of electron loops of the corresponding diagrams 
contributing to the photon vacuum polarization function. Thus we have 
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where N is (for formal purposes) the number of types of lepton. 
Substituting Eqs. (5) into Eq. (4) and Eq. (3), we get 
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Using the concept of the QED invariant charge 
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one can find similar relations between the 3-loop expressions for the photon vacuum 



1 The /5-functions and coefficients of Eqs. (1) and (2) are normalized as in Refs. [^, [22 
27| , not as in Refs. fl, |8|, |16| . 



polarization function in the OS scheme, namely U(x,a), and that in the MS scheme, 
renormalized at // = m, namely U(x,a). Substituting the large- a; expansions 

III = ( a i + b\ In x) , 111 = (oi + b\ In x) 

n 2 = (o 2 + &2 In x) , n 2 = (a 2 + b 2 In x) 

n 3 = (o 3 + 63 In x + c 3 In 2 x) 

n 3 = (a 3 + 63 In a; + c 3 In 2 x) (8) 

in Eq. (7), and taking into account Eq. (4), we get 



This means in turn that 
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(10) 
These identities are supported by the results of explicit calculations. (Compare the OS 



scheme results, summarized in Ref. 16 , with those of MS, obtained in Ref. M.) For the 



2-loop and 3-loop contributions to the constant terms of II(x, a) and II(x, a) we have 
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In Eqs. (10) and (11) we used notations for the coefficients of the photon vacuum polar- 
ization function similar to those defined in Eq. (5). Substituting Eqs. (11) into Eqs. (6), 
we get the following relations: 

p? ] = ?r i -M 1| (4 1| -4 11 ) 
flf = 3f-2/3f(4"-4 11 ) 

Pf = ?f-2/3f(4 2| -4 21 ), (12) 



which can be also obtained from the analysis of Ref. JT6| . 

3. Let us now consider the results of our calculations. The analytical expressions 
for the 2-loop and 3-loop coefficients of the photon vacuum polarization function in the 
MS scheme, II(x,a), were presented in Ref. |§. We are interested in the following gauge- 
invariant results [§]: 
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They were obtained in the Feynman gauge in the process of carrying out the work de- 
scribed in Ref. |T§ with the help of the SCHOONSCHIP program MINCER ]T2J. In the 
course of our work these results were confirmed by independent calculations made in an 
arbitrary covariant gauge with the help of the REDUCE program SLICER [B3|, especially 
written for this purpose. 

The coefficients g 2 , g 3 , and g 3 , were calculated in Ref. f| from 2- loop and 3- 
loop massive propagator-like integrals with zero external momentum, with the help of the 
REDUCE program RECURSOR. They read 
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These numbers were also confirmed in an arbitrary covariant gauge, by means of the 
specially written FORM program SHELL3 ||20| . 

Using the results of Eqs. (13) and (14), we get the following expressions for the 
constant contributions to the photon vacuum polarization function in the OS scheme: 
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The result for a 2 is in agreement with that known for many years |2T| . The expression 
for a 3 coincides with that recently obtained in Ref. [22], which is known to be consistent 
with numerical calculations [Q of 4-loop corrections to the muon anomalous magnetic 

m E21. 



moment a. 



The result for a 3 is new. It is obtained by combining results from Refs. |5], |8| that 



have been exhaustively checked by our two new programs |T9|, |2D] . Using SLICER |I!|, 
we calculated the bare photon self-energy, IIb, to three loops at large — q 2 . Then we used 
SHELL3 [|(| to calculate LTb to three loops at q 2 = 0. The constant Z 3 that gives the 
OS renormalized result 11(0, a) = also gave the asymptotic OS constants of Eq. (15) 
directly, without any reference to the MS scheme. 

Note the presence in a 3 of £(2) and C(2) In 2 terms, which are typical of calculations 
in the OS scheme. The origin of the ((2) In 2 term is well understood: it comes from 
a scalar massive integral that also contributes to the 2-loop correction to the fermion 
anomalous magnetic moment || |9], [QJ. Notice also that the £(5) contribution to the 
3-loop expression of U(x,a), like the 2-loop £(3) term, is independent of the choice of 
scheme. (Compare Eq. (15) with Eq. (13).) This fact can be understood by using the 
methods of Ref. [H . 
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strongly differs from the numerical result 
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obtained from the analysis of the results of numerical calculation of the 4-loop corrections 
to a M |TJJ by the methods of Ref . |23] , taking into account the correct analytical expression 



for the a 3 terrr^]. 
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4. By using our new analytical expression for the a 3 term, we can also obtain 
the analytical expression for the set of 4-loop diagrams contributing to a M formed by 
inserting the sum of 15 one-particle irreducible sixth order single-electron-loop vacuum- 
polarization diagrams in a second-order muon vertex. (See Fig. 1.) Indeed, using the 



methods of Ref. [23], which agree with those of Ref. ||24|| , we get 
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where Iq = 1/2, I\ = —5/4, and b 3 = 1/32 |25| . Substituting the analytical value of 03 
from Eq. (15), we have 
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whilst the numerical calculations of Ref. |14| gave 
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The discrepancy between Eqs. (19) and (20) is much greater than any revealed by the 

analytical calculations in Ref. |2(J of some 8th order contributions to the electron anomaly. 

The discrepancy between Eqs. (16) and (17) also leads to a substantial change to 

the 10th order contributions to a M of the diagrams of Fig. 2. Indeed, using the methods 



of Ref. [23], one can get the following expression for them [lG] 
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Taking the value of the a 3 term from Eq. (15), we are now able to obtain the following 
totally analytical expression for Eq. (21): 
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The corresponding numerical result reads 
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It should be compared with the similar expression obtained in Ref. |16| from Eq. (21), 
after taking into account the 'old' a 3 value from Eq. (17), namely 
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Therefore, one way of resolving the discrepancy between Eqs. (16) and (17) might be an 
accurate numerical calculation of the diagrams of Fig. 2. 

5. Our next step is the determination of the 4-loop correction to the QED (3- 
function in the OS scheme, in analytical form. We first remark that the results of the 
analytical calculations of the 4-loop corrections to the QED /5-function in the MS scheme, 
namely to the /3-function defined by Eq. (1), are 
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Using now Eqs. (5) and (6), and the results of Eq. (14), we get 

fa = -— A - - A 2 . 
H 16 9 



(25) 
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At iV = 1, Eq. (26) coincides with the result (3 3 = -121/144 obtained in Ref. J27]. 
The 4-loop coefficient f3± has the following form: 

p 4 = _^iv + fl_^ C (2) + ^C(2)ln2--C(3)W 2 
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where the powers of iV serve merely to distinguish contributions with corresponding num- 
bers of electron loops. (With leptons of different mass, the finite parts of OS charge 
renormalization involve dilogarithms ||, and the OS expansions of Eq. (8) involve log- 
arithms Jl7| of mass ratios; both are absent from the MS scheme.) Setting N — 1, we 

obtain 

5561 _ 23 8 7 

H 5184 9 SV ; 3 SV ; 8 SV ; ' V ; 

which may also be obtained from the m- independence of Aq| d , as given in Ref. 0. The 
order N term in Eq. (27) is scheme independent. The A 3 contribution coincides with 
the result of Ref. ||22|| . The A 2 contribution is new. In view of the results described in 
Sections 3 and 4 it differs from the one obtained previously taking into account the 
'old' numerical value of the a 3 term. (Compare Eq. (16) with Eq. (17).) 
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Let us now compare the numerical behaviour of the 4-loop approximations of the 
/3-functions in the MS and OS schemes with the corresponding function in the MOM 
scheme, namely with the \l/-function, which we here normalize in correspondence with 
Eqs. (1) and (2), as follows: 

*«mom = — ^j /3{a) = _ (3(a) , 29 

a In a a In a 

where «mom is the invariant charge of Eq. (7), at fixed q 2 = —A 2 . The expressions for the 
functions (3(a) and ^/(a) can be obtained from the results of Ref. 0: 

(3(a) = 0.667 g) + 0.5 g)' - 0.215 g)* - 1.047 g)' ... (30) 

a\ f a\ 2 fa\ 3 fa 



¥(a) = 0.667 (-j+ 0.5 y + 0.100 (-j - 1.202 (-j ... (31) 

In the OS scheme the analogous series reads 

(3(a) = 0.667 (-) + 0.5 (-Y - 0.840 (-Y - 1.142 (-V ... (32) 

Notice that now, in all three cases, the 4-loop coefficients are negative and have the 
same order of magnitude. This is a welcome result, in view of the natural expectation 
that the physical conclusions should be less scheme-dependent after taking into account 
higher-order perturbative corrections. 

Finally, we use the new numerical value of the 4-loop coefficient of the OS /9-function 
of Eq. (32), to obtain a slight modification of the estimate presented in Ref. |L6|| for the 



renormalizat ion-group constrained [E3| 10th order contributions to a^. Our result reads: 

o{, 10) (R.G.) = [B 5 + 82.92(78)] (-) (this work) , (33) 

where B 5 is the unknown 5-loop mass-independent contribution (expected to be positive) 
and the uncertainty in the numerical value of the renormalizat ion-group determined mass- 
dependent term derives largely from the uncertainty in B^ = —2.503(55) [[i||]. Our result 
differs little from that of Ref. [TJ, namely a{, 10 >(R.G.) = [B 5 + 86.57(78)] (a/rrf, obtained 
using Eq. (17). 

6. In conclusion: the 4-loop and 5-loop muon-anomaly contributions of Eqs. (19), 
(23) and (33) follow from our new results for the 3-loop OS vacuum-polarization coefficient 
ajj of Eq. (15) and the 4-loop OS /5-function coefficient of Eq. (27), each given in analytical 
form, for the first time. In the course of deriving a 3 = 0.3268745 ... we have devised two 
entirely new 3-loop programs, SLICER |19[ and SHELL3 [BO], which substantially validate 



the programs MINCER JL2] and RECURSOR [g], respectively, and hence confirm the large 
discrepancy between our result and the numerical estimate a 3 ~ 1.356 that follows from 



the 8th order muon-anomaly contribution of Eq. (20), obtained by Kinoshita et al. [I 1 



In this connection, we remark that there exists an impressive body of work |3Qj to support 
our belief that the methods of dimensional regularization, used here and in Refs. ||, |8|, [12[] , 
yield final results free of any infra-red pathology. Our value for 03 does not affect the 



result of Ref. [16| for the on-shell constant 64, whose N = 1 value has recently been 



verified [ESI. A formula involving 03 has been given for 65 PJ|. This is in error by the 



omission of a term involving the 5-loop coefficient in the expansion of Eq. (31), whose 
value is unknown. 



Note added: After completing this work, we communicated our result of Eq. (19) 
to Professor Kinoshita, who undertook a re-evaluation of the diagrams of Fig. 1, ob- 
taining |3]J a numerical result substantially different from that given in Eq. (20) and 



considerably closer to ours. 
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